Transport properties of overheated electrons trapped on a Helium surface 
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An ultra-strong photovoltaic effect has recently been reported for electrons trapped on a liquid 
Helium surface under a microwave excitation tuned at intersubband resonance [D. Konstantinov et. 
al. : J. Phys. Soc. Jpn. 81, 093601 (2012) ]. In this article, we analyze theoretically the redistri- 
bution of the electron density induced by an overheating of the surface electrons under irradiation, 
and obtain quantitative predictions for the photocurrent dependence on the effective electron tem- 
perature and confinement voltages. We show that the photo-current can change sign as a function 
of the parameters of the electrostatic confinement potential on the surface, while the photocurrent 
measurements reported so far have been performed only at a fixed confinement potential. The exper- 
imental observation of this sign reversal could provide a reliable estimation of the electron effective 
temperature in this new out of equilibrium state. Finally, we have also considered the effect of the 
temperature on the outcome of capacitive transport measurement techniques. These investigations 
led us to develop, numerical and analytical methods for solving the Poisson-Boltzmann equation in 
the limit of very low temperatures which could be useful for other systems. 



PACS numbers: 



i.20.Hh, 89.75.Hc, 05.40.Fb 



Electrons trapped on the liquid helium surface were the 
first experimental realization of a high mobility two di- 
mensional conductor P, Q • The extremely high mobility 
of the surface electrons allows to explore in an unique way 
many problems in fundamental physics, examples include 
Wigner crystallization [3|, propagation of magnetoplas- 
mon waves 0, [3| , transport of correlated charges in con- 
fined geometries [g, 7], quantum melting [8j] and sliding 
0. The interest in this system has been renewed by its 
potential for quantum computation that comes from their 
very large spin coherence times 10] and from the spacial 
addressability of the electrons on the surface Re- 



cently, an unexpected transport regime was observed in 
novel experiments where surface electrons were excited 
by a Millimeter- wave irradiation 13|, This regime 



occurs when the perpendicular magnetic field lies within 
regularly spaced intervals for which the ratio between 
the microwave and cyclotron frequencies is close to an 
integer value. It also appears only when the photon en- 
ergy is equal to the energy spacing between the two low- 
est transport subbands. This transition is usually called 
inter-subband resonance, it corresponds to transitions be- 
tween the ground state and the first excited state of the 
electronic wavefunction in the direction perpendicular to 
the Helium surface. Once the above conditions on the 
perpendicular magnetic field and on the photon-energy 
are satisfied, the irradiation can lead to a complete sup- 
pression of the dissipative component of the electronic 
response in capacitive measurements. This effect is strik- 
ingly similar to zero-resistance states that were first re- 
ported in ultra high mobility GaAlAs hetero-structures 
15, 16]. 

Further experiments on surface electrons showed a 



strong redistribution of the electronic density compared 
to its equilibrium shape under zero-resistance conditions 
[l7| . This redistribution was characterized by a deple- 
tion of the electronic density at the center of the electron 
cloud and by an accumulation of electrons towards the 
edge of the cloud. These experiments were motivated 
by the possibility of trapping at the edges of the system 
under microwave irradiation that was proposed theoret- 
ically [l8| . However, experiment and theory are still in 
their early stages and other mechanisms could be respon- 
sible for this redistribution. In order to distinguish be- 
tween the different mechanisms that can compete with 
each other we have analyzed theoretically the change in 
the electron density induced by a strong overheating of 
the surface electrons. 

When surface electrons (SE) absorb microwave irradia- 
tion their effective temperature can increase significantly 
above the temperature of the Helium bath. At zero mag- 
netic field, the effective temperature of the electrons un- 
der microwave pumping of the intersubband resonance 
was studied both theoretically and experimentally in 
It was shown that electrons could be overheated to tem- 
peratures in the range of Tg ~ 10 K even if the temper- 
ature of the Helium bath was much lower Tq ~ 0.3 K, 
however no estimation of Tg is available in the zero resis- 
tance state (ZRS) regime. Experimentally it was shown 
that the formation of ZRS is accompanied by a strong re- 
distribution of the electrons on the Helium surface, which 
can cost up to a 100 meV charging energy per electron 
17|. Thus the electrons need to absorb many microwave 
photons (typical frequencies are in the 100 GHz range) 
to store the needed energy, which suggests that they can 
occupy highly excited states. It is probable that their 
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energy distribution cannot be described appropriately by 
an effective temperature in this regime. However, since 
this is the simplest starting model to understand the dis- 
tribution of the electrons under irradiation, we have in- 
vestigated the effect of a high electron temperature on 
the electron density profile ne(r) in the trapping poten- 
tial. As the charging energy cost is surprisingly high we 
are led to consider temperatures much larger than the 
thermal energy of the bath, we note however that these 
temperatures are still an order of magnitude lower than 
the penetration barrier into the liquid Helium or than 
the height of the confining potentials which are both in 
the eV range. 



I. SIMULATION OF HOT ELECTRON DENSITY 
DISTRIBUTION 

Since SE form a non degenerate system, their distribu- 
tion on the surface is governed by a Boltzmann statistics 
at the electron temperature Tg: 
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Here e is the absolute value of the electron charge, ^(r) 
the electrostatic potential created by the confining elec- 
trodes and the electron cloud. The quantity is the 
total number of electrons trapped in the cloud. The 
electrostatic potential V{t) is determined by solving the 
Poisson equation in a cylindrical cell that is sketched on 
Fig. [TJ The potential and electron density depend only 
on the distance r to the cell axis and the potential V{y) 
can be expressed in the following form: 

V{r) - VeAr) + I G{r,r')n,{r')2nr'dr' (2) 



where Vext{r) is the potential created by the confining 
electrodes alone, in absence of the electron cloud, and 
G{r,r') a Green- function giving the potential created by 
electrons located at radius r' . Analytic formulas for 
yext{f) and the Green-function are derived in the ap- 
pendix for the case where the effect of the finite radius 
of the experimental cell can be neglected. This approx- 
imation is accurate in the limit where the difference be- 
tween the radius of the experimental cell and the radius 
of the electron cloud is larger than the cell height 2h 
(see Fig. [1]). Finite elements simulations in a more re- 
alistic geometry confirmed that finite cell size effects are 
expected to be negligible in the experimental geometry 
used in H, 14, 17|. The appendix also gives an explicit 
analytic formula for the potential Vext{f) as a function of 
the potential of the bottom disc B and guard electrodes 
Gi,2- The combination of Eqs. ()ll2p . forms a Poisson- 
Boltzmann equation that needs to be solved in order to 
determine the electron density Ueir) and the electrostatic 
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FIG. 1: Schematic diagram of the experimental cell modeled 
in our simulations. The cell is cylindrically symmetric, sur- 
face electrons are trapped above the central disc electrode B 
and the bottom guard electrode Gi. The potential of B was 
fixed to Vb = 4.2 V in all the simulations (a typical experi- 
mental value) . Electrodes C2 and Ci are grounded, except in 
section IV where a small AC potential is applied to Ci. In 
all simulations the potential difference between Gi and G2 is 
fixed to Vb for reasons that are explained in the text. The 
voltage Vgi is thus the main control parameter for the shape 
of the electron cloud. 



potential V{r). Interestingly we have previously studied 
a similar problem in the context of counterion condensa- 
tion around permeable hydrophobic globules [20l] . 

The numerical solution of this equation in the limit of 
low (but finite) temperatures is challenging since small 
errors on the value of the potential V{r) can lead to 
a large errors on the electron density ne(r) as a con- 
sequence of Eq. IT]). This can introduce instabilities 
in many numerical procedures (for example direct itera- 
tion of Eqs. pi2l) or relaxation methods based on drift- 
diffusion equations). From our numerical experiments, 
we found that the best stability at low temperatures was 
achieved by a Monte-Carlo computational method. 

In this method, an ensemble of Nt ~ 10'' trial particles 
were displaced on the simulated region of the Helium 
surface, using a Metropolis algorithm in the potential 
landscape given by Eq. The unknown density ne(r) 
that appears in this equation was determined by aver- 
aging over time the probability of presence of the trial 
particles in thin circulars shells in which the surface of 
the electron cloud was separated (typically the available 
configuration space was divided into 100 circular shells). 
The outlined numerical method was stable for tempera- 
tures as low as Te = 0.1 K, convergence was checked by 
computing the relative error: 
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with Z = Jd^r'exp ^^^^^ the normalization factor in 
the Boltzmann-distribution. In this expression the po- 
tential V{r) was computed by numerical integration of 
Eq. ^ for the electron density obtained from the Monte- 
Carlo algorithm. Note that a more detailed description of 
the numerical method that we used here will be published 
elsewhere 21| . The value of this relative error parameter 
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FIG. 2: Electronic Density profiles (r) at the liquid helium 
surface as a function of the radius r from the center of the cell, 
represented at two different electron temperatures (Te = 1 K 
and Te = 500 K) and for different potentials of the electrode 
Gi (Vgi = 3V : blue curves ; 5V : green curves and 5.4V : 
red curves) with a fixed total number of electrons A^e = lO'^. 
The potential of the lower central electrode B is Vb = 4.2 V. 



in our simulations was around eerr — 10 ^ at high tem- 
peratures (Te ~ 100 K), terr =i 5 X 10"^ at ~ 1 K and 
terr — 4 X 10~^ at the lowest temperatures Te ~ 0.1 K; 
as discussed above t^rr is higher at low temperatures due 
to the longer convergence times. 

Typical electron density distributions are shown in 
Fig. [2] for several values of the gate parameters for 
Nf, = 10^ electrons in the cloud. As the potential in 
the guard electrodes is increased we observe a transition 
between different shapes of the electron cloud. In the 
limit of a low potential on the guard electrodes, the elec- 
trons are confined in the central region of the cell below 
the electrode Ci and the shape of the electron cloud is 
described by a monotonic distribution where the density 
vanishes at large radius r. This case corresponds to the 
curves at Vgi = 3 V which are represented on Fig. [51 The 
value of the bias voltage between the guard Gi and the 
disc B is then Vqi — Vb = —1.2 V. As in the experiments, 
a potential is also applied to the top guard electrode G2 
in order to enforce the relation Vgi — Vg2 = Vb = 4.2 V. 
This constraint allows to keep the perpendicular compo- 
nent of the electric field constant across the cell which 
is important to keep the same intersubband transition 
energy for all the electrons in the system. On the con- 
trary, when the potential of the peripheral electrode Vgi 
is much higher than Vb, the electrons are confined only 
above the guard electrode Gi, forming a ring (see the 
curve at Vgi 5.4 V and 1 K on Fig. [2]). Finally an 
intermediate case appears when the Coulomb repulsion 
between electrons or their thermal energy is sufficiently 
strong to overcome the energy barrier created by the pos- 
itive bias potential Vgi — Vb- In this case the electrons 
are mainly confined under the electrode G2 but a part of 
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FIG. 3: Electronic density variation Afie ~ 

[n (Te) ~ n {Te — O.IK)] as a function of the guard po- 
tential Vgi for different effective temperatures (from top 
to bottom at Vgi = 1 Volt: O.IK, 50 K, 100 K, 200 K and 
300 K). The plain black curve corresponds to the predictions 
of the analytic theory Eq. Hllll2p at effective temperature 
Te — 300 K. In the inset, the density variation is rescaled 
by the temperature. For Vgi < Vg all curves collapse on 
each other with a very good agreement with the analytical 
predictions. For Vgi > Vc, the scaling breaks down since the 
photocurrent becomes thermally activated (see discussion in 
the text). 



the electrons density is also localized above the central 
electrode B forming a non-monotonic disk shaped dis- 
tribution ne{r). This case is represented on Fig. [5]) for 
Vgi = 5.0 V and Te = 1 K or 500 K, and Vgi = 5.4 V and 
Tf, = 500 K. Note the transition between a ring shaped 
density to a non-monotonic disk which is observed when 
temperature increases at Vgi — 5.4 V. 

The comparison between low and high temperature 
curves for these three cases reveals the following trend. In 
the case of a monotonic ring density the electron cloud 
tends to expand outwards as the electrons are heated, 
however when electrons form a non-monotonic disc den- 
sity the trend is reversed and the electrons tend to come 
back towards the center as Te increases. This effect can 
be qualitatively understood as the tendency of the elec- 
tron temperature to smooth the electron density profile 
in the cloud. 

In order to study quantitatively the effects of electron 
heating on the electron density we have focused on a 
quantity that can be accessed directly in experiments. 
In 



17| the change of the electron density was probed by 
measuring the transient currents created on the top elec- 
trodes Gi , G2 , Gi by a periodic change of the density pro- 
file induced by an On/Off modulation of the microwave 
power under ZRS conditions. Since the electrodes Gi 
and G2 are held at the same potential ne(r) is almost 
constant below these two electrodes, thus a simple plane 
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capacitor model can be used to relate the transient cur- 
rent ii{t) measured by a current amplifier connected to 
Ci and the change in the SE density: 
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where S{Ci) denotes the Helium surface below the elec- 
trode Ci. Within the frame of our simplified model, 
the On/Off modulation of the microwave power changes 
the electron temperature between its equilibrium value 
To ~ 100 mK and its (unknown but possibly several or- 
der of magnitude larger) steady state value Te under illu- 
mination. The time-integral of the photo-current during 
a microwave ON half-cycle will thus be given by: 



i,{t)dt = --S{Ci)ineiT,) 
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(5) 



where fie denotes the average electron density below the 
central electrode Ci. Thus knowing Ane = ne{Te) — 
ne(Tb), it is possible to estimate the amount of charge 
displaced in the experiment. The dependence of Afie 
on the guard voltage is displayed on Fig. [3] for several 
temperatures Tg. 

This figure shows that Afie changes sign as a function 
of the guard voltage Vgi, and exhibits a sharp maxi- 
mum for Vgi close to the transition from a ring to a 
non-homogeneous disc which occurs at Vqi = Vc with 
Vc — 5.21 Volt, this value is obtained from numerical 
simulations at low temperatures. For potential Vgi < 
Vb, a temperature increase causes an expansion of the 
cloud below the electrode Vq\ leading to An^, < 0. The 
quantity Afie changes sign at Vgi ~ Vb = 4.2 Volt (more 
precisely at a slightly higher value Vgi — 4.45 Volt) and 
when the bias voltage Vgi — Vb becomes positive we ob- 
serve Afie > 0. We note also that the density change Afie 
almost vanishes when SE form a ring with a vanishing 
density in the center at Vgi > Vc- The latter behavior 
can be explained with the following argument: at zero 
temperature when Vgi > Vc the potential at the center 
of the cell is lower than the potential inside the electron 
ring. This potential difference creates an energy barrier 
that precludes the electron from exploring the center of 
the cell when the thermal energy is small compared to 
the barrier height. The maximal values of Afie are ob- 
tained when the electron cloud is in a non-homogeneous 
disc configuration Vb < Vgi < Vc- This is intuitively 
plausible since in this case the center of the cell and the 
ring are at the same potential at Tg = and there is 
no energy barrier to prevent hot electrons from exploring 
the center of the cloud. As expected Afie increases with 
temperature. 

The above heuristic arguments can explain the sign of 
iWLe as a function of Vgi , but fail to account for some sur- 
prising aspects of the functional dependence such as the 
seemingly non-analytic behavior at Vgi — Vc- Indeed 



Afie (Vgi) features a sharp maximum for Vgi < Vg, but 
quickly vanishes at Vgi > Vg, the asymmetry of this 
peak becomes even more pronounced at lower Te- In 
the next section, we develop a perturbation theory that 
allows to explain this singularity and its temperature de- 
pendence. 



II. PERTURBATIVE CALCULATION OF THE 
EFFECT OF TEMPERATURE 

In this section, we derive a simplified analytic theory 
to explain the behavior of the SE when the electronic 
temperature is raised from to Te. At zero temperature, 
the potential of the electron cloud is fixed and equal to 
a constant value Vq. This allows to effectively linearize 
the equations since the electron cloud can be replaced by 
an equipotential surface. The boundaries of the cloud are 
then fixed by the constrain that the electric field vanishes 
on the boundary (stability condition) and the potential 
of the cloud Vq is fixed by the number of electrons in the 
cloud Ne - Thus the problem of finding the zero tempera- 
ture density distribution ne{r,T — 0) is computationally 
much simpler, and we assume that ne{r, T = 0) is known 
in the following derivation. 

If the temperature changes, the number of electron in- 
side the cloud changes by an amount of Sue, modifying 
the potential inside the cloud by an amount 5V - We 
assume that the density inside the cloud is just a little 
perturbed Sue ^ ne{r,Te = 0): 

ne{r,Te) = ne{r,Te^Q) + 5ne{r)- (6) 

Using the Boltzmann-equation Eq. ([1} we also have: 
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where Uav and n^^, are normalization constants. Using 
our assumption Sue ^ ne (r, Te = 0) we are led to 



ne{r,Te = 0) ~ riat, exp 



/ eSV 



(8) 



\kBT, 
or equivalently: 

eSV = kBTe [log Tie (r, T = 0) - log Uav] - (9) 

The quantity 6V can be estimated in an other way; in 
the center of the electron cloud where the density is uni- 
form, the variation of density Sug leads to a change in 
the potential SV: 
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this expression just comes from the electrostatics of a 
planar capacitor. The value of Sue taken at the center 
of the cell is equal to Aue, eliminating 5V between the 
equations © and ([TU]), we obtain the following relation 

^_ 2eo/ci3Te[logne(r = 0,r = 0)-lognJ , , 
Ane = . (11) 

This equation is accurate only in the center of the elec- 
tron cloud where density gradients can be neglected. 
Outside the electron cloud, the density decays exponen- 
tially with the temperature and can be neglected. If we 
tentatively assume that Eq. (jlip is valid in the entire 
cloud, we can find Uav from the requirement of the con- 
servation of electron number: J Srierdr = 0, leading to 



/cwlogK (r,r = 0)]rdr ^ 

J Cloud 



(12) 



The integral in this expression is performed for the elec- 
tron density at T = for which the limits of the cloud 
are well defined. 

The combination of Eqs. (Illll2p allows to find the vari- 
ation of the electron density at the center of the cell based 
on the knowledge of the electron density at T = 0, they 
also imply that Ane scales linearly with the temperature 
Te (we assume Tg is much larger than the Helium bath 
temperature To). The above analysis is justified only for 
< Vc where the density He {r — 0,T — 0) is non van- 
ishing, we have thus compared our semi-analytical the- 
ory with the simulations only in this domain. The results 
are displayed on Fig. [3l and show that our semianalyt- 
ical approach is very successful for Vg — Vc where it 
reproduces the numerical results with very good accu- 
racy without any adjustable parameters. The agreement 
is worse in the limit Vq ^ Vb probably because our as- 
sumption that Eq. ()lip is valid in the entire cloud does 
not hold anymore when the electrons are strongly con- 
fined by the guard electrodes. The scaling Afie oc Tg is 
also confirmed by the simulations as it is illustrated in the 
inset of Fig. [S] To summarize the proposed theory seems 
to explains quantitatively the change of the electron den- 
sity under heating, and allows detailed comparison with 
numerical simulations and hopefully future experiments. 



III. EFFECT OF HEATING ON TRANSPORT 
MEASUREMENTS 

In previous sections we have investigated the effect of 
an elevated out of equilibrium temperature on the dis- 
tribution of the electrons on the Helium surface. The 
information on electron density can be accessed through 
photo-current measurements similar to those performed 
in [l3|. However the most wide-spread type of experi- 
ments on SE is a Sommer- Tanner technique which allows 
to determine the electron longitudinal mobility fi^x p^ . 



We have thus analyzed numerically the effect of Tg on 
the outcome of such experiments. 

In a Sommer- Tanner measurement an A.C. potential 
Vac cos uit is applied to one of the top electrodes (we have 
taken this electrode to be Ci in our simulations) an the 
induced AC displacement of image-charges is detected 
from the remaining ones using a lock-in technique. Our 
numerical approach to simulate this experiment was a 
direct numerical integration of the time dependent drift- 
diffusion equations: 



due 



f^xx 9 

r dr 
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(13) 



This equation is coupled to a Laplace equation on the 
electrostatic potential, which includes the electron den- 
sity and the static and time dependent potentials on the 
electrodes as a boundary conditions. The drift-diffusion 
equations were integrated using an implicit first order 
integration scheme solved with a finite element method 
(FEM). For reference we give here the weak formulation 
of the differential equation connecting ne{r,t + A) and 
ne{r,t) in our numerical calculations: 

[ne(r, t + At) — Heir, t)] u{r)rdr 



f f , , ,dV du 



dV du ksT due (r, t + At) du 



dr 



dr 



rdr 
(14) 



the potentials are evaluated at time t and recomputed at 
time t + At {At is our integration time step) using the 
updated value of the density rie, u{r) denotes a finite ele- 
ment trial function. Finite element simulations were per- 
formed using the FreeFem programming language |23| . 
for typical parameters — 3 x 10^, Vb — 4.2 V, 
Vgi — Vg2 = V. The AC potential with amplitude 
Vac — 3 mV was applied to the electrode Ci, and the 
modulation of the image-charge number was determined 
on the neighboring electrode C2 after subtraction of the 
direct capacitive coupling. This modulation Q{t) was 
separated into in-phase and out of phase components: 
Qx and Qy which respectively denote the amplitude of 
the cos and sinwt components in the total signal Q{t). 
An immediate consequence of Eq. (1131) is that the quan- 
tities Qx and Qy depend only on the ratio this is a 
consequence of the fact that the potentials react instan- 
taneously to the displacement of charges at the typical 
experimental frequencies in the kHz range. This quan- 
tity should then be compared with the parameters of the 
electron cloud which, for a monotonic disc density dis- 
tribution are the density at the cloud center and the 
cloud radius Rc- These considerations lead us to define 
the dimensionless parameter W = '^'^'^■^•^ which charac- 
terizes the ratio between the frequency and the response 
time of SE. 

We summarized our numerical results in a Cole-Cole 
diagram that displays trajectories followed by Qx and Qy 
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FIG. 4: Cole-cole diagram of the in-phase (Qi) and out of 
phase (Qy) charge modulation on electrode Ci as the modu- 
lation frequency oj is changed. The AC potential is applied 
to Ci and the charge modulation due to the direct capacitive 
coupling between both electrodes was subtracted. The contin- 
uous curves were obtained at different temperatures that are 
indicated in the legend (as temperature increases the semicir- 
cle radius decreases), the direction of increasing lo (decreasing 
/in) is indicated by the arrows on the curves. The dashed 
curve with symbols joins points at different temperatures for 
a fixed value of u)/ = 2 cm^/V. 

on (Qx, Qy) plane as the measurement frequency uj is var- 
ied at fixed Te . At low frequencies when W the elec- 
trons have time to follow the time varying potential adi- 
abatically and all the displaced image-charge is in phase 
with the excitation Qx ^ Qy As frequency increases, an 
out of phase component appears since the electrons are 
unable to respond fast enough to follow the phase of the 
external potential. Finally at frequencies much higher 
that the electron characteristic response time ^ 1 
their response vanishes and the trajectory converges to 
the origin {Q^ = 0, = 0). 

The simulation results at different temperatures are 
shown on Fig. 21 the obtained trajectories in the {Qx,Qy) 
plane are very similar to semicircles. As the temperature 
is increased the radius of the semicircle tends to shrink, 
qualitatively this occurs because some of the terminally 
exited electrons are not sensitive to the applied bias po- 
tential. As a consequence the shape of the trajectories 
in the {Qx,Qy) plane, are very different if the mobil- 
ity of the electrons changes under microwave irradiation 
or their temperature. Indeed when fi^x is changed the 
points {Qx, Qy) still collapse on their semi-circular trajec- 
tory because it is parametrized only by the ratio . On 
the contrary, when the temperature changes the points 
{Qx,Qy) lie on different semicircles which leads to very 
different trajectories. An example of a trajectory versus 
temperature at fixed oj/^xx = 2 cm^/V is shown Fig. SI 
Thus heating can create deviations of the {Qx,Qy) tra- 
jectories from the behavior that would be expected if only 
the electron mobility changed under microwave irradia- 



tion. 



IV. CONCLUSIONS 

We investigated the effect of an elevated out of equilib- 
rium electron temperature on the distribution of surface 
electrons on liquid Helium in a typical experimental ge- 
ometry. Our main result is probably the determination of 
the photocurrent induced by an increase of the electron 
temperature under irradiation. This was achieved nu- 
merically by developing a Monte-Carlo method capable 
of computing the electron density even in the limit of very 
low temperatures, and analytically through a perturba- 
tion theory analysis of the effects of temperature. The 
dependence of the photocurrent on the control gate volt- 
ages that is predicted within our model is quite specific 
and should provide a reliable signature of electron heat- 
ing, against other mechanisms possibly at play. Namely 
we predict that the photocurrent generated by density 
variation at the center of the cloud changes sign as a 
function of the guard voltage and is stronger in the re- 
verse bias configuration. Unexpectedly, we also find that 
the photo-current is maximal when center of the electron 
cloud is almost depleted. Finally we have also demon- 
strated numerically that in Sommer- Tanner experiments, 
heating should induce trajectories on the Cole-Cole di- 
agram that differ strongly from those expected from a 
change in electron mobility. This could be checked in 
transient conductance measurements, where the trajec- 
tory on the Cole-Cole plane can be measured as the sys- 
tem is driven in and out of a zero resistance state by an 
on/off modulation of the microwaves power. We hope 
that our results will allow a quantitative analysis of the 
electron energy distribution under zero-resistance con- 
ditions, and that the analytical and numerical methods 
developed here will be useful in other systems that are 
described by a Poisson-Boltzmann equation. 

We are thankful to E. Trizac, D. Konstantinov and K. 
Kono for fruitful discussions. One of us, A.C. acknowl- 
edges the support from RIKEN and St Catharine college 
in Cambridge. 



APPENDIX 

In this appendix, the electrode potential Vext (r) and 
the Green function G (r, r') are calculated. The total 
potential V (r) is the sum of the electrode potential 
Vext (r) and the contribution coming from the electrons 
= I G (r, r') Ue (r) d'^r. The cylindrical symmetry 
of the system allows us to write the different potential 
as functions of the cylindrical coordinate r = {r,9,z) but 
we also make use the Eulerian coordinate r = {x,y,z). 
Note that in our notations ||r|j ^ r. 
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Electrode potential Vext (r) 

The potential Electrode Vext (r) is solution of the Pois- 
son Equation 



0. 



(15) 



Taking the xy- Fourier transform of this equation ()15|) . 
we get : 



dzzVext - k^Vext = 

where k'^ = + ky, k — {kx, ky) and 
d^k 



Vextir)= II —e^'^-'^VextiKz) 



(16) 



(17) 



The solution of Eq. pB)) are 

Vext {kx,ky, z) = A{kx,ky) cosh(fcz) 
+B{kx, ky) sinh(fcz) 



(18) 



where A (kx, ky) and B (kx, ky) are functions which sat- 
isfy the bordering conditions. The cylindrical symmetry 
imposes A (kx,ky) — A (k) and B (kx, ky) — B (k), lead- 
ing to 

/>oo 

Vext (r, 0,z) ^ / fcdfc Jo (kr) [A(k) cosh(fcz) 

S(fc) sinh(fcz)] , (19) 



J 



where Jj (x) is the i — th Bessel function of the first 
kind. Quoting Vi the potential of the electrode i, Ri 
the maximal radius of the electrode i and H (r < Ri) the 
Heaviside function (equal 1 if the argument is true, else 
zero) , the electrode potential has to respect the following 
boundary conditions : 

Vext(r, z^h) = Vg, [H(r < Rg,) - H(r < i?c.)] (20) 
Vext(r, z^-h) = VeHir < Rb) + 

Vg, [H{r < Rg,) ~ H{r <Rb)]. (21) 

It is straightforward to calculate 

/•OO 

H(r < R)= RJi(kR) Jq (kr) dk. (22) 
Jo 



From the equations (l20l) . ([21]) and ([22|) . we obtain the 
following boundary conditions : 

Vext{r, h) 

= /o" Mkr) Vg, [RG.MkRG.) - Rc,Ji{kRc,)] dfc 
= /o°° kJn{kr) [A(k) cosh[fc/i] -)- B(k) sinh[fc/i]] Ak 
Vextir, —h) 

= ir Mkr) [(Vb - Vg,) RBJi(kRB) 

+VG,RG,Ji{kRG,)]dk 
= /o°° kJoikr) [A(k) cosh[kh] - B(k) smh[kh]] dk . 

Which leads to: 



A(k) = 



1 



[Vg, [RG.MkRG,) - RcMkRc,)] + (Vb - Vg,) RsMkRB) + VG,RG,Ji(kRG,)] , (23) 



2k cosh[kh] 

^ TnrAjTu^^GARG,Ji(kRG,)-Rc,Ji(kRc,)]~(VB~VG,)RBMkRB)'VG,RG,MkRG,)]- (24) 
2k smh[fcn,J 



r 



With these expressions of A(k) and B(k), the elec- 
trode potential Vext (r) is completely determined. 



Electron potential Ve (r) 

Lets us consider an electron in the center of the sys- 
tem, that means localized at the coordinate (0, 0, 0). This 
electron creates a potential Ve (r). This potential can be 
written 



^e(r) 



47ren 



V 



(25) 



The term l/||r|| is coming from the solution of the 
Laplace equation for a punctual charge and the term 
V, solution of the Laplace equation AV = 0, per- 



mits at the potential to satisfy the boundary conditions 
Ve (r, 6, z ^ ±h) = 0. As previously, the potential V can 
be obtained as, see Eq. p^ : 



V(r,z) 



kdkJo (kr) [A(k) cosh(fcz) 

-f i?(fc) sinh(fcz)] . (26) 



The function B (k) must be equal to due to the sym- 
metry of the problem. Furthermore, we can calculate 



1 f 

]M\ ^ fcJo(M— ^dfc, 



(27) 
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leading to the full electronic potential Ve (r) : 



47reo Jo 



kjf){kr) 



A{k) cosli[fcz] 



(28) 



The boundary conditions Ve {r, z = ±/i) = lead to the 
following expression for A (k) : 



A{k) 



k cosh[kh] 



and then 



47reo Jo 



Joikp) 



coshffczl 



cosh [fc ft] 



(29) 



dfc. (30) 



The complete potential V (r) 

To obtain the complete potential due to all electrons 
and electrodes, we have to sum the electrode potential 
and the potential created by the surface electrons dis- 
tributed with a density rie (r) in the plan z = 0: 

V[r,z) = Vext{r,z) 

Vedlr - r'll , z) ne(r') r'Ar'AB. (31) 



Now, we have to note that the electrons can move only 
in the plan z = and, by this fact, the electrons feel a 
potential y(r, z = 0) = V(r) defined by : 

V^(r) -l/e,t(r,z = 0) 



- — / dfctanh[fcft] / drV'Jo(/cr) J^ikr') n(r') . 
2eo Jo Jo 

(32) 

The comparison between Eq.(l32|) and Eq.© give the fol- 
lowing expression for the Green function 



G(r,r') = 



47reo Jo 



d/ctanh[fcft] Jo(fcr) Jo(fcr'). (33) 



The full potential V{r), see eq. ([32]) . is numerically ob- 
tained. In order to increase the numerical precision, the 
integrations are done successively between two consecu- 
tive zeros of the Bessel functions Jo{kr) or Jo(fcr'). 
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